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^ I The use of parametersHrneasuring order-parameter fluctuations (OPF) has been encouraged by the 

C . recent results reported inp which show that two of these parameters, G and Gc, take universal values 

fy>j ' in the liniT^o- In this paper we present a detailed study of parameters measuring OPF for two mean- 

field models with and without time-reversal symmetry which exhibit different patterns of replica 
symmetry breaking below the transition: the Sherrington-Kirkpatrick model with and without a 
field and the Ising p-spin glass (p=3). We give numerical results and analyze the consequences 
which replica equivalence imposes on these models in the infinite volume. We give evidence for the 
transition in each system and discuss the character of finite-size effects. Furthermore, a comparative 
1-^ I study between this new family of parameters and the usual Binder cumulant analysis shows what 

(~| , kind of new information can be extracted from the finite T behavior of these quantities. The two 

Q ' main outcomes of this work are: 1) Parameters measuring OPF give better estimates than the 

O , Binder cumulant for Tc and even for very small systems they give evidence for the transition. 2) 

For systems with no time-reversal symmetry, parameters defined in terms of connected quantities 
are the proper ones to look at. 
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One of the main features of the low-temperature phase of spin-glasses is the lack of self-averageness of the order 
parameter. Parameters measuring order-parameter fluctuations (OPF), A, A^ G, and Gc, where introduced to locate 
phase transitions in systems which did not contain rtime-reversal symmetry (TRS) in the HamiltonianEl. Very recently 
the use of these parameters has been enhanced ina (henceforth referred to as RS) where it has been shown that G 
'"^ ' and Gc take, in the low-temperature limit, the universal values 1/3 and 13/31 respectively, provided the ground state 
P5 , is unique and that there is no gap at zero field. This situation is generically met in spin-glasses with no gap at zero 
O ■ coupling, thus bringing more interest to the applicability of these parameters in the study of phase transitions. 

In ordered systems, a good parameter to locate phase transitions is the Binder cumulant which is the kurtosis of 
the order-parameter probability distribution. The uniqueness of the ground state in that case is enough to guarantee 
that the Binder cumulant takes the universal value 1 at zero temperature for any finite volume. On the contrary, in 
, ■ , order to find universal values for Gc and G, one has to add a condition, namely the absence of gap in the local field 
j^ ' distribution. 

Concerning disordered systems, it is well known that the use of the Binder cumulant may present serious problems, 
specially in models without TRS in their Hamiltonian. Then, the Binder cumulant may take negative values showing 
no evidence for a crossing of the different curves for different sizes. Recently, Hukushima and Kawamurafl have shown 
that this is specially true in models with one-step replica symmetry breaking where the Binder cumulant may take a 
negative value just below Tc. This explains the non-positiveness of the Binder cumulant in Potts glasses and p-spin 
models. The same being also true for p even when there is time-reversal symmetryQ. 

However, at finite T, parameters which measure OPF always remain finite and positive. In systems with non- 
vanishing OPF, the behavior of these parameters in the thermodynamic limit can be inferred from the restrictions 
that the replica equivalence property imposes on the P{q), q being the overlap between two different replicas of the 
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same system. Even though the behavior of G is trivial, the behavior of Gc is not, and depends on which kind of 
repUca symmetry breaking pattern (RSB) characterizes the frozen phase. This knowledge can help us to distinguish 
among different systems and can also provide us with useful information about the spin-glass (SG) phase. 

One of the main outcomes of RS is that already very small sizes may yield valuable information about the infi- 
nite vQlauae spin-glass transition. This ig confirmed by previous xiumerical studies of Ising spin glasses and p-spiii 
modelsld'EHS, Migdal-Kadanoff spin glassesEl, Heisenberg spin glasses!. Potts glassesEl and some protein folding modclsEl. 
In view of the main results obtained in RS we will concentrate our attention on models where the new parameters 
G and Gc-, as well as A and Ac may yield valuable information not attainable using the Binder cumulant analysis. 
In particular, we are interested in models without time-reversal symmetry such as spin glasses in magnetic field and 
p-spin models. 

Our purpose is to present a detailed numerical analysis of parameters which measure OFF by studying very small 
samples. Numerical calculations have been performed using exact computations of the partition function for small 
samples and Monte Carlo simulations. Monte Carlo calculations have been done in the spirit of the study on the 
spherical Sherrington-Kirkpatrick (SK) model shown in RS where very small sizes were simulated. One of the main 
outcomes of RS is that calculations of G and Gc at low-temperatures need a very large number of samples. This is 
due to the fact that at low temperatures only those samples having local fields much smaller than T give a significant 
contribution to G and Gc, so that the uncertainty in determining the value of G and Gc becomes very large. In other 
words, the low-temperature behavior of G and Gc provides direct information about the quality of the data. Having 
values for G greater than 1/3 shows either a lack of thermalization or poor sample statistics, a situation which the 
Binder cumulant does not reflect, since it is much less sensible to these fluctuations. 

The paper is organized as follows. The next section presents the definitions used in this paper. Section III analyzes 
the consequences the replica equivalence imposes on the values of these parameters at finite T. Sections IV and V 
analyze two mean- field models exhibiting different kinds of RSB: the SK model with and without a magnetic field 
and the Ising p-spin model with p=3. Finally, in section VI, we present the conclusions. Two appendices are devoted 
to several technical points. 

II. A SHORT REMINDER 

Order-parameter fluctuations are an intrinsic characteristic of disordered systems. Although extensive quantities 
are self-averaging the order parameter may not. This is the typical situation found in systems where due to RSB 
there is a multiplicity of states and ergodicity is broken. Also in the case where replica-symmetry is not broken OFF 
are very sensitive to the spectrum of fluctuations around the stable solution. Spin glasses and disordered systems in 
general are characterized by the presence of large scale domain excitations without a typical length scale, and as a 
result correlations have a very slow decay, typically as a power law, and hence, with no characteristic length scale. 
The equilibrium phase is then marginally stable and OFF decay slower than 1/V where V is the volume. The possible 
scenarios concerning the behavior of OFF have been discussed in RS. Here we content ourselves with defining the 
parameters of interest in the analysis of this paper. 

Four parameters which measure OFF were studied in RS: G, Gc, A, Ac, where the subindex c stands for connected 
disorder quantities (to distinguish from disconnected quantities). B, the usual Binder cumulant and Be, its connected 

insional spin glass with binary couplings. Calling q — y X^i SfS^ 
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counterpart were used ir£j to study the four dimensional spin glass with binary couplings. Calling q — y X]i '^i'^-, 



the spin overlap of two systems of spins 5f and S^, with V the number of spins, we define these parameters as follows, 

(1) 




A=^'^'^'_JP , (2) 



G=<f^-^,,^'-^, (3) 

(9^) - {q") 



{{q-{q)f) ) 






{{q-{q)W-{{q-{qm 1 A, 

{{q-{qm-{{q-{qm ^ 1 - Be 

All of these parameters vanish at high temperature. It is wellinown that in the limit T — > the Binder ratio takes 
the universal value 1, provided the ground state is unique. Iro it was shown that under the hypothesis of a unique 
ground state and no gap in the local field distribution these parameters satisfy the following properties, 

^imG(T/,T) = i (7) 

]imG.{V,T)^^ (8) 

lim A{V, T) = . (9) 

The last identity for A follows trivially from (g) due to the uniqueness of the ground state. In Appendix A we make 
a brief summary of the main outcomes of section VI in RS to show how these universal values are obtained. The 
low-temperature results for Be and Ac are also derived in that appendix. The results are: 

lim AciV,T) -> oo , lim Bc{V,T) -^ -oo . (10) 

These last two relations follow from the fact that even though numerator and denominator for both quantities vanish 
at zero temperature, the denominator vanishes as a higher power of T, and thus the ratio diverges. 

We have to note that A and Ac have a non trivial behavior so they must be analyzed in detail in order not to extract 
wrong conclusions at low temperatuces as shown by B. Drossel et al. in a very recent study of these parameters in 
the Migdal-Kadanoff approximationtll. In particular, one has to be careful when dealing with transitions at T = 0, 
as the limits T ^> and y ^ oo do not commute, so that results in (|o|) could not hold. l-. 

All these quantities are dimensionless, thus we expect that a finite-size scaling analysis can be performedQ: 

X = fx (l{T ~ TcY'^'^) X = G,A,B. (11) 

Quantities which can be obtained from derivatives of the free energy, i.e. thermodynamic quantities, are expected to 
have the same critical exponents, hence, for the Binder cumulant we have that vb — i^ (the usual correlation length 
exponent). However, we do not know a priori how to express A^G in terms of thermodynamic quantities, so we 
cannot assure that scaling exponents for G will be the same as the ones for B. The equivalent of (0) definedin terms 
of the magnetisation instead of q,-has been extensively studied by numerical simulations for diluted modelalj and for 
the Ashkin- Teller random modelli3. In both cases the exponents associated with the parameters A and B {G is only 
function of A,B, see eq. (0)), within numerical precision, have been shown to be the same. Here we reach similar 
conclusions for the SK model. 

Despite that critical behavior for A, B, G is governed by the same exponents {i.e. by the critical behavior of the 
same correlation length) the parameters which measure OPF appear to be good indicators for a phase transition. 
In the spin-glass phase G and Gc attain a finite value and vanish in the paramagnetic phase. On the contrary, A 
may vanish in both the spin-glass phase and paramagnetic phase but stays finite at Tc (a result predicted for diluted 
systemsEil). The interest of studying A is that it does not vanish when OPF are finite, a situation characteristic of 
replica symmetry breaking transitions. Such information can not be deduced only from G, Gc, Be and Ac which may 
be finite even when OPF vanishuH. 

Let us stress that only in cases where time-reversal symmetry is absent both class of parameters (connected or 
disconnected) can be studied. In case of models with time-reversal symmetry in the Hamiltonian only disconnected 
averages (G, A, B) are feasible. A calculation of the disconnected parameters GctAc and Be would require to 
artificially break time-reversal symmetry which may impose additional difficulties to the simulation. This explains 
why ina only the parameters G and A were studied for the Ising spin glass chain and the spherical SK model. 



III. REPLICA EQUIVALENCE 

Systems exhibiting RSB at the transition are well-known to have a non trivial P{q) below the transition temperature. 
The information about this probability distribution is contained-jp the replica matrix Qab, where sub-indices stand 
for different replicas. The replica equivalence property (RE)EjtZl states that every row/column of this matrix is a 
permutation of any other row. Thus, any quantity such as ^^^ Qab does not depend on a or, equivalently, the free 
energy must contain terms proportional to n (which staiids-|for the number of replicas) in order to ensure that there 
will be no divergences when performing the limit n -^ OlJ'Ej. 

This property imposes some constraints on the P{q) = Pj{q). i-Qne can derive general relations for the 
probability distributions of several overlaps and, as is pointed out iriij, any probability function of m overlaps 
P(gi2, (?34, (756, ..., 92(m-i)2m) Can bc cxprcsscd in terms of the probability distribution of one overlap P{qi2) plus 
the cyclic distributions, i.e. pi2, 23,31^^12,23,34, 41^^^ pi2, 23,34,. ..,mi^ For the case of two overlaps, fluctuations satisfy 
the following relation: 

1 2 

P{qi2, 934) = ^P{qi2)5{qi2 - qzi) + ^P{qi2)P{qzi) , (12) 



which in the thermodynamic limit yields the well-known Guerra relation: 
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{<ih)'-^{qt2) + ^{Qh) , (13) 



I.e. 



which, as stressed in RS yield for systems with non- vanishing OFF (n3) a trivial behavior of G, 

G = ie(T - T,) . (14) 

In RS it was pointed out that this result could hold even in spin-glass systems with a marginally stable replica 
symmetric phase. Notwithstanding, there is no such general trivial behavior for the other disconnected quantities, A 
and B, neither for Gc, Ac or Be. To compute Gc one has to deal with the probability of functions depending on three 
and four overlaps, so that in general one has a much more complicated object which depends on the following terms: 

^^Y.Q^.-Jn,),^dq, (15) 

a,b 
—f Tw Ty: y] QabQbcQca = / ^^(^12, 923, q3l)qi2Q23q3ldqi2dq23dq31 (16) 

n[n- l)(n-2) ^ J 

TiT, 57 V" QabQbcQcdQda= / ^(<7l2 , g23'?34, g4l)'?i2'?23'?34'3'4lrf<?12C?g23rf934rf<741 



n(n — l)(n — 2)(n 

^ a,b,c,d 



TrQ" 



i(n-l)(n-2)(n-3) 



(17) 



Therefore, we expect that at finite T the behavior of Gc will depend on the kind of RSB that takes place at the 
transition, as it depends strongly on the structure of Qab- When computing the connected quantities, Ac and Be 
one encounters the same situation so that we have to go to particular patterns of RSB in order to obtain a simpler 
expression for Gc, Ac and Be, as well as for A and B. 

In particular, for the case of one-step RSB, where Qab can take two different values qo and gi, and the n replicas are 
divided in m blocks, we obtain a very simple expression for all the connected quantities (calculations which lead to 
these results are shown in Appendix B): 

_ 39- 113TO + 98m2 
" ~ 93-22lTO + 266m2 ' ^^ ^ 

_ 39-113m + 98TO^ 
14077i(l — m) 

3(31- 167m +182m^) 
^^ 280m(l-m) " ^^^^ 



These results establish that connected quantities depend only on the RSB parameter m, which in general is a function 
of T. However, A and B depend also on the elements of the replica matrix, namely go and qi, as shown in the 
expressions ( [79| ) and ( |80| ) in Appendix B. 

The simple functional dependence on m of Ac, Gc and Be lets us establish several universal features of these 
parameters. In first place, we note that below Tc, Gc always remains finite and positive and that at T = 0, since 
TO = in any system, it takes the universal value 13/31d. In the same way, at to, = 1, Gc takes another universal 
value: Gc = 0.174. Another interesting feature of Gc is the existence of a local minimum at to, = 0.650 whose value 
Gc = 0.113 is universal for all the systems exhibiting one-step RSB. As we will show in a forthcoming section, this is 
important from the point of view of numerical results for small sizes as the asymptotic position of this minimum can 
give us information about the behavior in the thermodynamic limit and thus about the dependence of m on T. 

In second place, we note that both Ac and Be diverge at to, = 0, and therefore at zero temperature, the same 
being also true at to = 1. The difference relies on the fact that while Ac remains always positive and has a minimum 
at m — 0.56 whose universal value is Ac = 0.187, Be diverges to — oo and displays a maximum at a positive value 
Be — 0.316 ( TO, — 0.451). This is an interesting result, on the one hand because Gc, Ac and Be provide us with 
valuable information about the characteristics of the SG phase, and on the other, because the exclusive dependence 
of the connected quantities on rn establishes a direct connection with models for structural glasses, where m, the RSB 
parameter, is well-known to control the violation of the fluctuation-dissipation theorem (FDT). 

Unfortunately, for the two-step RSB scheme (and further RSB), we do not have such a simple expression and 
numerical computations for each system have to be performed. In the two following sections we analyze two mean- 
field models: the SK model in a field and the p-spin {p = 3) Ising model. Neither of the two models contains TRS 
in the Hamiltonian and exhibit RSB of different kinds: a full step and one-step RSB respectively. Therefore, these 
systems are characterized by having non- vanishing OFF and can provide a direct check of the results obtained in this 
section. 

IV. THE SHERRINGTON-KIRKPATRICK MODEL WITH AND WITHOUT A MAGNETIC FIELD 

This section studies the Sherrington-Kirkpatrick model with and without a field. Preliminary results for the model 
in a field wjere shown inH. For the model without a field results have been recently presented by Hukushima and 
Kawamurau. Here we present a more systematic and detailed study for small sizes of connected and disconnected 
quantities. The SK model in a field is defined by the disordered mean-field Hamiltonian, 

HsK — - 2_^ Jij '^i ^3 ^ ^ z_^ ^j • (21) 

i<j i 

Before discussing the results in a field we discuss the case of zero field. An important diffcrsace between the study 
of OFF parameters and Binder cumulants (for very precise recent SK simulation results seeEilu3) is that very high- 
precision statistics is needed to compute OFF parameters. Consider for instance G and Gc- These are ratios of two 
quantities which may be very small if OFF vanish yielding a large error for G and Gc- In models where OFF vanish 
simulating large sizes may then require a prohibitive computational effort. Not only long simulations are needed 
to thermalize the samples but high-precision statistics is needed to determine with reasonable precision the ratio of 
quantities which vanish in the infinite- volume limit. This second limitation is not present in models where OFF are 
finite such as the SK model at finite temperature. 

A. Zero-field simulations 

Results for G and A are shown in figures |l| and g. In figure |l|, we display results ior^^ve^iy small sizes {N — A — 11) 
obtained from exact computations of the partition function (such analysis was done inE3 for the SK model) averaging 
over 10000 samples of Gaussian quenched couplings with varian ce || l /A^. In figure g we show results for larger sizes 
{N — 32, 64, 128, 256, 512) using the parallel tempering techniqueESEj with binary couplings. The number of samples 
ranges from 1000 for the smallest size, N = 32, to 250 for the biggest one, N = 512. Figure g shows the Binder ratio 
B for both groups of sizes. In figure we show the results of the scaling analysis performed to the results for B (a) 
and G (b) for the larger sizes. |-| 

Monte Carlo simulations for larger sizes corroborate the results for G and A reported inQ. It is important to note 
that, as this system contains TRS, thc—Biiider cumulant for different sizes is expected to exhibit a crossing of the 
curves for different sizes as reported irJjO. However we must note that, while even for very small samples we do 
observe a crossing for G and A at T close to Tc = 1 (figure |l|), for the same sizes, curves for B do not display such 



a crossing (figure || (a)). In figure ^ we have plotted the Binder ratio for A^ = 32 together with the smaller sizes to 
stress that we have to increase the size of the system {i.e. reach N « 100) to observe the crossing and not only a 
merging of all the curves in the low T region. Results for larger systems display a crossing of the curves of all the 
parameters as expected. Our observations corroborate previous resultsQ which argued that finite-size corrections to 
Tc for parameters measuring OPF were of opposite sign than for the Binder parameter. For G and A (figures m and 
H) the crossing point at T* starts well above Tc and approaches Tc as we increase the size of the system. For the 
biggest sizes T* is very close to Tc. On the contrary, curves for B cross at T* < Tc, and therefore, since thermalization 
is more difficult at low T, the crossing is harder to observe. Moreover, the crossing for the biggest sizes takes place 
at a temperature slightly below Tc, thus making it evident that finite-size corrections to Tc are much stronger in the 
Binder ratio than in A and G. 

The finite-size scaling analysis shown in figure shows a remarkable resuljU, In the particular case of the SK model, 
the argument of the scaling function ( |ll| ) for the Binder ratio can be shownE3 to be A^ (T — Tc)^. The collapse of the 
data, which are shown only for T — Tc > 1, is very good, not only for B but also for G (the same being true for the 
parameter A, although not reported here). This suggests that exponents for G are the same_as for B, justifying a 
search for a general field-theoretical argument in the line of that proposed for diluted systemalJ. 

As far as the details of the calculations are concerned, we have to point out that for small systems Gaussian 
couplings assure the uniqueness of the ground state, so that as reported in RS we obtain the universal values for G 
and Gc. For binary couplings, however, this results should not hold for small systems, but as long as we increase 
the size of the system (i.e. we approach the thermodynamic limit) wc should not expect a dependence on the bond 
distribution because in mean-field entropy vanishes at T = also in models with discrete couplings, and therefore, 
universal results should also hold. In order to be sure that results can be compared, we have checked that already for 
A'^ = 33, 32 results for both coupling distributions nearly coincide. 
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FIG. 1. Parameters A (a) and G (b) for the SK model with no field by exact computations of the partition function averaged 
over 10000 samples for A^ = 4, 5, 6, 7, 8, 9, 10, 11 from bottom to top at low T. Error bars are shown for sizes A'^ = 4, 7, 11. In 
both insets we show in detail the crossing region for sizes N = 4,7, 11. The solid line in (a) corresponds to numerical results 
of the one-steo-approximation (|79|), obtained from solving numerically the saddle-point equations of the SK model for this 
particular casecj. 




FIG. 2. Parameters A (a) and G (b) for the SK model with no field for the case of binary couplings for sizes 
iV = 32, 64, 128, 512 from bottom to top in the high T region. 
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FIG. 3. Binder ratio, B, for the SK model with h — 0: (a) for small systems {N — 4, 5, 6, 7, 8, 9, 10, 11 from top to bottom) 
and (b) for iV = 32, 64, 128, 512 from top to bottom at high T. 
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FIG. 4. Scaling functions of B (a) and G (b) for the SK model with ft = for sizes iV = 32, 64, 128, 512. 



B. Finite-field simulations 

In a field one must be a bit careful. According too two conditions are necessary to ensure equations (Q) and (||): 
uniqueness of the ground state and no gap in the local fields distribution. For the SK model with Gaussian couplings 
these conditions are immediately satisfied. But the first condition may be violated in the SK in a field. This happens 
in samples which have a ground state with zero magnetization. In this case the ground state is twofold degenerated 
because both the zero magnetization configuration and its fully reversed one have the same energy. In computing 
OFF for finite volume these rare samples have a finite probability violating the requirement of uniqueness of the 
ground state. Consequently eq. (pt) for Gc is not satisfied. Note that, on the contrary, eq. (R) is still satisfied because 
it is invariant against time-reversal symmetry. A zero magnetization ground state can only occur for samples with an 
even number of spins, so we expect strong differences for cases where N is even or odd for small systems. 

We have performed simulations with two applied fields h = 0.3 and,ii.= 0.6. The transition into the frozen phase of 
the SK model in a field is determined by the de Almeida-Thouless linecSl which yields T^^q^ — 0.65 and T^^^g = 0.48. 
In figure g, we show results at field h — 0.3 for an odd number of spins for sizes TV = 3, 5, 7, 9, 11 by exact calculation 
of the partition function and N = 13, 17, 23, 33, 64 by Monte Carlo simulations. The largest size {N = 64) has 
been simulated with binary couplings and parallel tempering, while the smallest ones have been simulated by simple 
Metropolis algorithm with Gaussian distributed Jij . As we have said at the beginning of this section, when systems 
are small and the ground state is not unique, eq. (^) is violated as can be observed in the inset of the plot for Gc 
(figure 1^ (c)) where Gc becomes greater than 13/31, and goes to 1 at T = 0. This does not happen when systems 
are small and there is an odd number of spins. Results for Gc and G for the largest field, h = 0.6 are shown in 
figure ^ for sizes A^ = 7, 11, 13, 23, 33. Results for Ac, the Binder ratio and its corresponding connected quantity. Be, 
are displayed in figures M, M and O, for both applied fields, h = 0.3 (a) and h = 0.6 (b). 

• Small fields: h = 0.3. 

First of all, let us focus our attention on the results for h = 0.3. In figures H and R we observe that G, Gc, A 
and Ac display a clear crossing of the curves which is not seen neither for the Binder ratio where curves for all 
sizes seem to merge in the low temperature (figure pi (a)), nor for the connected Binder cumulant (figure ffa)). 
Despite, there is a clear difference between disconnected parameters (G and A) and connected ones (Gc and 
Ac). The crossing for disconnected quantities, A and G, takes place at a temperature which is higher than 1 for 
small samples and slowly approaches 1 as we increase the size. On the contrary, results for Ac and Gc, display 
a crossing, T* , at a temperature which already for small samples is smaller than one, and that approaches the 
predicted Tc = 0.65 as we increase the size of the system. Indeed, for the largest sizes A^ = 29,33,64 curves 
cross around T — 0.65 (figure |^ (c)). 




FIG. 5. Results for the G,A, and Gc for the SK model in a field h — 0.3. a) Curves for G for A'' = 5,6,7,9,11- dashed 
lines -, 17, 23, 32, 64 from bottom to top at low T; the full line corresponds to the infinite-volume limit result G = 1/3. b) Results 
for A, for sizes: N = 5,7, 9, 11, 13, 23, 33 (from bottom to top looking at the curves around T ~ 0.75). c) Parameter Gc for 
sizes: A^ — 3,5,7,9,11,23,32,64. (Note that in c) there's no full line standing for the T = value of Gc as it falls out of 
scale.) In the inset on the top right corner, we show the low T behavior of Gc for small samples with an even number of spins 
N = 2,4, 6, 8, 10 from top to bottom, the full line corresponds to the T — value Gc = 13/31. Curves below the line correspond 
to samples with odd number of spins A'^ = 3, 5. 



• Large fields: h = 0.6. 

Results in a larger field, h — 0.6, for G, Gc and Ac displayed in figures ^ and (b) show a greater difference 
between G and Gc- While the crossing for G takes place at very high temperatures (close to 2) , curves for Gc 
cross at a temperature which even for small samples is close to Tc = 0.48 and gets closer to it as we increase the 
size. From the Binder cumulant (figure p^(b)) one gets no clear information. Numerical results display a crossing 
at high temperature as curves for G do, but the crossing point seems to be moving towards higher temperatures, 
which suggests that this crossing has no real connection with the existence of a transition. Still, curves for Be 
(figure ^ (b)) show an interesting result, since we observe a crossing around T* — 0.5. Unfortunately, results 
are not very clean because the transition takes place at a very low temperature and thermalization is not easy. 
For this reason the crossing of curves for Ac in figure u\ (b) is blurijed by the existence of the divergence at 
T = 0. Notwithstanding this, our results, together with the results ino for the P{q), which for very large sizes 
presented two peaks in the positions predicted by the RSB solution, arc the most clear evidence for a transition 
in such a high field {h = 0.6) which has been reported up to now. 




FIG. 6. Gc (a) and G (b) for the SK model in a field h = 0.6 for sizes iV = 5, 7, 9, 11, 17, 29 from bottom to top in the low T 
region. 
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FIG. 7. Ac for the SK model in applied field for sizes N = 7, 11, 13, 23, 29: (a) for h = 0.3, (b) for h = 0.6. 
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FIG. 8. Binder ratio, B, for the SK model in applied field for sizes iV = 7, 11, 13, 23, 33, (a) h = 0.3 and (b) h = 0.6. 
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FIG. 9. Binder connected ratio. Be, for the SK model in applied field for sizes N = 7, 11, 13, 23, 33, with h — 0.3 on the left 
and /i = 0.6 on the right. 



C. Summary of the results 

• SK at zero field: A and G are good parameters to locate the transition. Finite-size effects are evidenced by 
a crossing at a temperature sliglitly above Tc, that moves towards Tc as size increases. The location of the 
transition point is more precisely located by G and A than by the Binder ratio, B. Finite-size scaling analysis 
reveals that B, G and A have the same scaling exponents, thus suggesting that the scaling behavior of A and 
G for generic short-range systems may be obtained from the standard renormalization group approach. 
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• SK in a field: In general connected quantities, Gc, Ac and Be show a crossing which gives evidence for a 
transition. Finite-size corrections to T^ are of the same sign as those in the case /i = 0. Fipte-size effects are 
stronger at low fields because small systems are strongly affected by the h — Q fixed pointed, suggesting that 
transitions in a field need to be studied at high magnetic fields. 

• Continuous RSB: In the frozen phase and for the largest sizes, the general temperature behavior of the different 
quantities, A, B, and Gc, Ac and Be in models in a field, is the same with and without a field. Therefore, it 
corroborates our results of the previous section: systems exhibiting the same kind of RSB below the transition 
(as in these two systems -the SK model with and without a field- which both display full-step RSB), display 
qualitatively the same behavior for A, B, Gc, Ac and Be- Note that we have intentionally excluded G from the 
list since it always takes the value (for V —> cxd) 1/3 in the low-temperature phase, hence does not depend on 
the particular RSB pattern. 

V. ISING P-SPIN, P=3 

In this section we turn our attention to a model without time-reversal symmetry where OFF are finite. We have 
studied the mean-field p-spin (p=3) model with Ising spins defined by: 

H = - 2_, Jii,i2,i3'^il'^i2'^i3 y (22) 

where couplings are Gaussiari_yariables with zero mean and variance 2^ ■ This model is well-known to have a one-step 
RSB transition at Tc ~ 0.6 Eil showing a discontinuity in the order parameter, which at the transition point jumps 
from zero to a finite value. These kind of transitions are driven by a collapse of the configurational entropy (also 
called entropy crisis), and have-been the object of a deep study as they have many points in common with the glass 
transition in structural glassesS. Particularly, we are interested in the fact that a measure of the violation of FDT 
in glasses can give us information on the pattern of RSB, i.e. on the parameter m which determines the structure of 
the replica matrix Qab- 

The interest of this model relies not only on the fact that it does not contain TRS in the Hamiltonian, but also on 
the fact that, as far as OFF do not vanish, results obtained in section III for the one-step RSB hold for this system 
and, thus can provide us with a direct check of the validity of identities (|l^), (|l^) and (gO|) for the finite T behavior 
of Gc, Ac and Be respectively. Moreover, this model presents a further simplification since qo, the off-diagonal block 
value of the replica matrix clement Qab, vanishes below the transition. For this reason, expressions ( |79| ) and (KG) 
given in Appendix B for A and B have a much simpler expression only in terms of m : 

A^—^, (23) 

3 (1 — m) 

2-3to 

B = — . 24 

2 (m - 1) ^' ' 

In order to study the behavior of OFF in this model, we have made exact numerical computations of the partition 
function for sizes ranging from A^ = 4 up to 11 averaging over 10000 samples and Monte Carlo simulations using the 
parallel tempering technique for N — 13, 17, 23, 29, 33 averaging over 150 — 500 samples. Results for B, Gc, G and A 



are displayed in figure nU and for Ac and Be in figure 11, the full line corresponding in each case to the theoretical 



prediction obtained by numerically solving the saddle point equations. We have to point out that despite the low 
number of samples due to the hard task of thermalizing these systems, results are strikingly good. 

Our results for the Binder parameter (figure Ru) do not display any crossing but show that around the transition 
temperature, Tc — 0.6, B becomes negative, and has a minimum which grows together with the size of the system 
as is expected from expression (EJ) which indicates that at the transition, whexe m is set to one, the Binder ratio 
has a negative divergence. Results obtained by Hukushima and Kawamura irjj for the Potts model, are different 
{B is expected to take the value —1 at the transition), since the actual values of qi, go and m are different, thus 
corroborating our result for B obtained in Appendix B (see eq. ^^, which states that B has a strong dependence on 
the model. Results for Be are very alike, since there is no crossing, and there is also a minimum around Tc which 
deepens with size. This behavior is very well understood by looking at the infinite volume behavior of both parameters 
in figures Iffl (b) and O (a) in which curves for different sizes approach to the theoretical curve from top to bottom 
and this is why there is no possible crossing. 
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Instead, our results for Gc^ G, A and Ac exhibit a crossing which moves from higher to lower temperatures giving 
evidence for the existence of the transition. As we had noted in the previous section, we observe that finite-size 
corrections to T* are bigger in disconnected averages than in connected ones, the formers giving a better estimate of 
Tc- However, it is important to note that the position of the maximum of A which accounts for the existence of a 
phase where replica symmetry is broken has a quick saturation towards T — 0.6, and grows with size as is expected 
from the divergence that this parameter shows at the transition (|2^). A similar situation is met in results for Ac 
(figure O), since there is also a predicted divergence at to = 1, which is numerically observable by a maximum around 
Tc which grows and sharpens with size. 

From results for Gc, we can still go a bit further. We observe that, as the size increases, the shape of the curve 
approaches the one of the thermodynamic limit much faster than in the other parameters, due to the softer behavior 
of Gc- It is interesting to note that even though we see a crossing of the curves which moves toward Tc from higher 
temperature as the number of spins increases, much more clear evidence for the existence of the transition and the 
kind of transition {i.e. one-step RSB) come from the existence of a bump which becomes closer to Tc as the size grows. 
We observe that at low-temperatures curves always remain under the infinite- volume line, so that there is no crossing 
in the low T phase. Moreover, curves quickly stuck to the value of the minimum displayed by the thermodynamic 
value of Gc. As we noted in the previous section, the fact that this value is universal can yield valuable information 
about the low-temperature phase of systems exhibiting one-step RSB. 

Results in the former section showed that connected quantities, Gc, Ac and Be for the SK model (figures || ^, |^, 
y) exhibited quite a different behavior, as the curves were smooth and did not show any concavity, nor divergence at 
the transition temperature. Hence, in addition to showing what kind of information we can extract from the behavior 
of Gc, Ac, Be, A and B, we also give clear ouantitative evidence of how different RSB transitions look like extending 
the previous comparative analysis made in u between the 3-Potts model and the SK with no field. 




1.5 



"^t^^« «i«»lM i] 



J , 1.5 

FIG. 10. We show the results for the 3-spin model for sizes ranging from TV = 4 to A^ = 33. In every plot, the solid line 
corresponds to the theoretical prediction of expressions (|l4), (|l8[), w%, and (pOl) for G, Gc, A and Be respectively (see 
text). On the left side we display, in the top plot, results for Gc for sizes N = 4,7, 11, 17,29,33 from top to bottom in the 
low T region, errors bars are showed for the larger and smaller sizes. Below we show the results for B (Binder ratio) for sizes 
N = 7, 11, 17, 23, 29, 33 from bottom to top . On the right side, we show results for the non-connected averages, A for sizes 
A^ = 5, 8, 11, 17, 29, 33 and G for sizes A'' = 4, 8, 11, 17, 29 both from bottom to top in the low-temperature region. 
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FIG. 11. Results for the 3-spin model: (a) Be for sizes A'' = 4,7,11,17,23,29 from top to bottom, and (b) Ac for sizes 
A^ = 4, 7, 11, 17, 23, 29 from bottom to top (in the low T region). The sohd hne corresponds to the numerically solved theoretical 
prediction for Ac and Be, see eqs. (uM and (GOl) respectively. 



VI. CONCLUSIONS 

In this paper we have presented a detailed analysis of parameters measuring OPF, A, Ac,G and Gc, as well as B and 
Be, for small systems. Despite the numerical effort that computation of these parameters may imply, we have shown 
the powerfulness of these parameters and what kind of information we can extract from them. Our analysis through 
replica equivalence in section III, shows that the intimate connection of OPF with RSB yields a dependence of A, B, 
Gc, Ac and Be on the RSB pattern, in opposition to G, which in the SG phase takes the value 1/3 at T < Tc ( even 
in marginal phases where OPF vanish, as stated in RS). Thus, A, B, Gc, Ac distinguish between systems displaying 
different kinds of RSB in its frozen phase. Moreover, we have established a direct connection with glasses, since we 
have shown that at any temperature below Tc, connected quantities depend exclusively on m (which determines the 
structure of the replica matrix Qat), which is the parameter having the physical meaning of an effective temperature 
which controls the violation of FDT in structural glasses. 

We have analyzed mean-field models which exhibit time-reversal symmetry in their Hamiltonian and models which 
do not have this symmetry. The universality of G and Gc at zero temperature stated in RS has been checked and has 
been extended to the rest of quantities under study A, Ac and Be and we have provided with a meaningful example of 
when identities ([7|) and (^ do not hold, i.e. the SK model in a field for small samples with an even number of spins. 
We have given many evidences that these parameters are useful tools to locate phase transitions into the spin-glass 
phase in any case. We can enumerate the main conclusions, 

1. Advantage of A, G (OPF) compared to B (non OPF). For models where an analysis of the Binder cumulant is 
enough to show the existence of a transition, i. e. the SK model without a field addressed in section IV, we have 
shown that already for very small samples, G and A give evidence for the existence of this transition. In other 
words, parameters measuring OPF have smaller finite-size corrections, so that we can extract useful information 
with no need of going to large systems. Moreover we have shown that A and G verify scaling relations and that 
the exponents coincide with those of the Binder ratio. Obviously, A, G require more computational effort than 
B but this effort compensates for the smallness of the corrections. 

2. Advantage of connected parameters compared to disconnected ones. For models with TRS only A, G, B are 
feasible to be computed. The interest of parameters involving connected averages concerns models which do 
not display TRS such as the SK model in a field and the Ising 3-spin model. For those cases we have once 
more checked the usefulness of these parameters to show the existence of a transition. A comparative analysis 
of G and Gc has shown that finite-size corrections are smaller for the latter, thus suggesting that connected 
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parameters such as Gc, Ac and Be are the best to look at in situations where there is no time-reversal symmetry. 
Our results for the SK model in a field show thatr-tjp, effects of the h = fixed point on small sizes are very 
strong (in agreement with the conventional wisdomc^O) , indicating the necessity to go to low temperatures- 
high fields to show the existence of a spin-glass transition in a field. In view of some recent numerical results 
this goal is probably not out of reach. 

3. Temperature dependence of OPF parameters on the RSB pattern. By studying these mean-field models which 
have different kinds of RSB transitions, we have checked that in systems where OPF do not vanish (i.e. when 
the results obtained through replica equivalence analysis are meaningful), Gc, Ac, Be, A and G have a different 
behavior in the low-temperature phase. In addition, we have shown that in the 3-spin model, which has 
a discontinuous transition, already for very small samples, Gc yields information about the infinite-volume 
behavior obtained in section III (eq. R9t) and thus about the replica symmetry breaking pattern through the 
temperature dependence of the parameter m. 

To sum up, our results for G, Gc, A and Ac for small systems, in addition to what had already been stated in 
RS, i.e. the universal values at T = for both G and Gc, and the trivial behavior of G for any system having a SG 
transition, establish that parameters measuring OPF meet all the requirements that a proper parameter to locate 
phase transitions should have in disordered systems. We have studied in detail their behavior in systems where OPF 
do not vanish and have shown what kind of information small samples give about the infinite-volume behavior at 
finite temperature. In conclusion, the study of parameters measuring OPF opens a new direction in the study of 
the SG phase, as they are useful not only to locate the transition, but also to establish the nature of spin-glasses 
at finite temperature. Moreover, the scaling analysis reveals that there might be a connection between OPF and 
critical fiuctuations that still lacks theoretical insight. The understanding of the trivial behavior of G in the infinite- 
volume limit can give us more information about the mechanisms that control this frozen SG phase and can help us 
to understand the connection between spin-glasses having a marginally stable replica symmetric phase and systems 
having a broken SG phase. The further extension of these ideas to non- disordered glass models is also an interesting 
open problem. 
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Appendix A: Low-temperature behavior of A, G, B and Ac, Gc, Be- 

In this Appendix we give a brief summary of the main points of the proof of the conjecture presented in RS, and show 
how the zero temperature values of all the parameters introduced in section II are obtained. We consider the most 
general finite spin system of Ising variables, described by the local field Haniiltonian: 



H = - ^ /izCTz , (25) 



under two hypothesis: 

• it has a unique ground state (up to a global flip of all the spins), given by the configuration {a*} : Eq = 

• its local field probability distribution has a finite weight at zero field. 

Under these two assumptions it was shown in RS that the low-temperature behavior of all the terms which appear 
in the definition of G and Gc (and hence in all the other quantities defined in section II) are dominated by one spin 
excitationsE3. This is because one-spin excitations yield a linear contribution in T, while two spin excitations yield a 
contribution of order (T^). Thus, in any general finite system where all excitations are possible, one-spin excitations 
give the dominant contribution at low T and are responsible for the universal values G = ^ and Gc = ■^ 'At T — 0. 
In cases where the lowest excitation systematically involves clusters of spins a generalization of the present argument 
shows that the present results still holdcJ. 

The argument proceeds as follows. First of all we note that all elements appearing in the definition of G and A can be 
expressed in terms of the two and four point correlation functions. If we define Tij = {(Tiaj)'^ and Tijki = {aiCFjak<Ji)'^ , 
then we have: 
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(«') = -^ + -^E^»^' (26) 



(9')' = ^ + 4E^'^^ + wE^'+4 E T..^T,^ + ^, E ^'^■^^'' (27) 



^2\ _ 



(«") = -[74 (3^' - 2^ + (6^ - 8) ^ T,, + 5] T,,fe,) . (28) 



J#J {i,j,k,l) 



Terms appearing in the definitions of the connected quantities also depend on one and three points correlation 
functions. With the definitions Ti = {ai)"^ and Tijk = (aiOja-kf' , we obtain the following expressions : 

(9)' = -^ (E ^' + 3 E ^'^' + 4 ^ TfT, + 6 Y. TlT.n) , (29) 

{q'){q) = -^ ((3^^ - 2) ^ Tf + (ay - 2) ^ t^t^ + 3 ^ tj^^.^ + ^ t„- ^r,) , (30) 

{q'){qf = ^ (^ E ^^ + ^ E ^^^^ + E ^^T^k + 4 E ^^^^'^^^ + 2 E ^%- + E ^^.^-^^-^O • (31) 

One-spin excitations. 

We consider now all possible one-spin excitations. The energy cost of reversing the spin i will be: e^ = 2h*a* . From 
here one can compute the correlation functions previously introduced and obtain: 

/ \ *f^ 2exp(-/3eO ^ 

^""'^ =^^['- l + E..,.exp(-2/3e0j " ^^ ^^ ^ '^^^ ' ^^'^ 

{a,aj)^a*a*{l-2{x,+xj)) , (33) 

{ai(jj(7k) ~ (J*(7*(7l (1 - 2(a;j + Xj + Xk)) , (34) 

{ai(jj(Tk(Ti) « (T*(T* (jla*i (1 - 2(a;i + Xj + x^ + x;)) , (35) 

with the definition Xi = exp(— /3ei). Note that we have made the approximation Xi « 1, thus we can approximate 



the denominator in (32) by one. We have to stress that, as widely discussed in RS, this approximation implies that 
we are in the limit f3V << 1. This is a key point in our argument concerning the local field probability distribution. 
When performing the averages over the disorder one has to tackle with the following objects: 



x" = / dh, dhj P{K,h.j) e-2/5(™ '»*'^*+« ''>*) (36) 



where the probability distribution of two fields results from integrating out the rest of the fields: 

Pihlh*)^ J \[ dhlP{hl,....,h*y) . (37) 

We have to point out that two local fields may be correlated so that the probability distributiaa. of one field P{h*) = Pi 
depends on the site. However, the probability of having two equal local fields is very small E3, thus in ( |37| ) there are 
no contributions of the type 5{h* — h*). It then follows from expression (|3^) by a simple saddle-point calculation 
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that terms which depend on two sites have a contribution ex 0{T'^) while terms which depend only on one site have 
contributions ex 0{T) : 



rF„ 0,0 TPJO) , , 

^ mn m 

By looking at the above expressions for the correlation functions one observes that terms which depend on only one 
site are present only in one-spin excitations since, as shown in RS, correlation functions computed in the case of 
two-spin excitations depend on the products XiXj^ thus bringing sub dominant contributions. 

Then, after some algebra, the results obtained for the quantities appearing in the definitions of G, Gc-, A^ Ac and 
B are the following: 



-2 



i6rE,^»(0), 



{q')'^{q') = ^"^^^:''"' iV-ir + Oir'), (39) 

W^W)W - ^'^^o^;f^"^ + ^(^') ' (41) 

W^JW) = ^-^f^+0{T'), (42) 
1 496r^^P,(0) 



((9 - il))') = ^ ^""iot;;""^ + ^(^') • (43) 
Then we obtain for G and Gc for any finite volume, 

G^l + 0{T) G, = H + o{T) . (44) 

The situation is rather different for Ac and Be- From the expression of eq. (0) of Ac, we see that the numerator is 

2 

linear in T, while the denominator {{q — (q))^) is quadratic in T. Hence we have that in the low-temperature limit: 

lim Ar (X — -^ oo . (45) 



The same situation occurs for the computation of Be given by eq. {m where we have to compare ((q — (q))^) ex T 

r2 



with {{q- {q})^) oc T^. It then follows that: 



lim Br oc — > — oo . (46) 



Appendix B: Using replica equivalence to compute OPF. 

In this appendix we report explicit calculations of some results outlined in Section III. In particular, we describe in 
detail how to make use of the replica equivalence property which states that quantities such as '^^ Q'^i^ do not depend 
on the replica index b, to compute the different terms involved in G and Gc which lead to identities (0) and (111) 
respectively. 

In general, we have to deal with terms involving one or more overlaps which in terms of the replica matrix elements 
Qab read : 

W) = -r^, T.Qab^ f dqP{q)q^ ^ q^ , (47) 

2^a bed l2( -1) h -1 ^O.b^cd--^l2y-l,l2y f 

(q)^ = ' 7 "'\u — H^— 7 77y 7T^ = / dqidq2 dqyP{qi,q2....qy)qiq2 qy , (48) 

n[n-\)[n-2)...[n-[2y-l)) J 

where sums always run over different indexes. Note that we have defined qk in order to simplify the notation. As we 
have already said in section III this can also be expressed in terms of the P{q) for a different number of replicas as 
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we show in the last equahty. 
Computation of G 

The terms which appear in the definition of G are the following: 

^ - n(n- 1) E Q'fa = / dqP{qW ^ g2 , (49) 

ab 

W) = -7-^ E Q'& = / ^9^(9)9' ^ 94 , (50) 

^ ' afc -^ 

^= n(n-l)(n-2)(.-3) ^g/°-Q-- = /^.c^.iP(.,.0.^.? ■ (51) 

Thus we only have to deal with terms involving one and two overlaps. As is shown iiO and discussed in detail inEj, 
the RE property allows to express the probability distribution of two overlaps in terms of the probability distribution 
of one overlap: P{q)- Hence (g^)^ can be expressed in terms of (g^) {q^) and (g). 

• Quantities involving two overlaps. 
The general term reads : 

^^^ - .(.-l)(.-2)(n-3) ^g^ Q^.Q^. • (52) 

The calculation proceeds as follows: the sum appearing in (p3) can be re-expressed as: 

E Q'l^bQ'cd = Y.Q'bY.Q'-b-2Y. Qab' - 4 E Q^-bQL ■ (53) 

a,b,c,d ab ab ab a,b,c 

In order to compute the last term of the previous identity we use replica equivalence, i.e. : 

c ac 

Thus, for the last term in (pSh, we obtain: 

E QabQL -T^Qabi: QL = E Qab (^^ QL) , 

a.h.c ab c ab 



(55) 



so that substituting this last result into identity (|51| ) we obtain the following general relation: 

J-EUZIT ^ {n-i){n- 1) Eab Q'ab Eab Qib , 2 Egfc Qaf 

^^ '^^ ' {n- 2)(n - 3) n(n - 1) n(n - 1) [n - 2){n - 3) n{n ~ 1) 

(n — 4)(?i — 1) 2 

= 7 ov 5Tft9' + 7 ^TT ^9fc+i • (56) 

(n — 2)(?i — 3) (n — 2)(n — 3) 



^From this relation, the computation of (g^)^ and terms such as (q) (g^) and (g)^ which appear in the definition 
of Gc is straightforward. Thus for (g^)^ we obtain: 

(n-mn-l) fEabQlb\\ 2 j:^,Ql, 



{q'){q')= 



(n - 2){n - 3) V"(" - 1)/ [n - 2)(n - 3) n{n - 1) 

(n-4)(n-l) 2 2 

(n-2)(n-3)'^2 + („ _ 2)(n - 3) ''' ' ^^^^ 
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At this point we are able to compute the value of the numerator and denominator of G, which yield the well-known 
result: 

numerator = — -( / dqP{q)q'^ - ( / dqP{q)q^f) 

denominator — i dqP{q)q^ — ( / dqP{q)q^)^ 

.. „ numerator 1 

-^ hm G = — - , (58) 

n^o denominator 3 

provided OPF (/ dqP{q)q^ — (/ dqP{q)q^)^) do not vanish. For the parameter A this yields in the limit n ^ Q : 

Computation of Gc 

Now we have to deal with harder terms which appear in the calculation of Gc and involve the probability distribution 
of three and four overlaps. The numerator reads: (g^)^ — 2{q'^){q)^ + (q)"' — ((g^) — {q)'^Y ^^^ the denominator reads: 
W) + 6(g2)(g)2 - 3"(^ - 4(173) (g) - (J^ - l^f so that we have to compute: 

WW)^^-, Tw .w ^ ^w 777 7T E QlbQcdQef (60) 

n(n — l)(n — 2)(n — 3)(n — 4)(n — 5) ■'^-^ 

('?)'* = —7 TTT ;^T7 ^77 TTT TTT ^77 ^ /, QabQcdQefQgh ■ (61) 

ra(n — l)(n — 2)(n — 3)(n — 4)(n — 5)(n — 6)(n — 7) ■'^-^ 

As we have already pointed in section III, these terms will depend in general in terms containing P{q) and 
terms which depend on the probability distribution of 3 or 4 connected overlaps, terms like J2a b c QabQbcQca ^^'^ 
Sa bed QabQbcQ^aQca whcrc rcplica equivalence cannot be used since we have in general that: ^^ QabQbc depends 
on the replica indices a and c, so that for any particular structure of the replica matrix Qab will have different values. 

• Quantities involving three overlaps. 
Now we will derive the general relation for any term containing three overlaps such as (pO[). The general term 
is: 

E QabQLQTf^ T. QabQLJ2QTf-i E Q'bQLQTb-^ E QabQLQZ- m 

a,b,c,d,e,f a,b,c,d e,f a,b,c,d,e a,b,c,d,e 

The first term can be computed from identity (p6|), since we have: 

-. TTT ,w ^ ,,, TTT r, E QabQidT^Qef = , "^"7 ^^ ^M ){q^){q^) ■ (63) 

n(n — l)(n — 2)(n — 3)(n — 4)(n — 5) ^-^ ^-^ ■' (n — 4)(n — 5) 

The other two terms on the r.h.s of relation ( |62|) are identical by permuting sub-indices {d, b) and super-indices 
(l, k) so that we only have to compute one of them: J2a bade QatQcdQeb- ^^ can apply RE (^) to the sum of 
index e and obtain: 

E QabQidQTb^ E Qa.Qcd(^^^)- E QarQL-2 E QabQcdQTb- (64) 

a.b^c^d,e a.b.c^d a.b^c.d a.b.c^d 

Again we make use of identity (pq) to compute the two first terms on the r.h.s of the previous identity, so that 
we are only left with the last term: ^^ j^ ^ ^ QtbQidQ^cb- -^y ^PPlyi^^g the RE to the sum over d we finally get: 

/ , QabQcdQcb — /_^ QabQcb I ' ~ I ~ 2^ QabQcbQca ~ /_^ QabQcb ' (^5) 

a^b,c,d a,b^c \ / a,b^c a^b,c 
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So, now, we know how to compute all the terms in function of the probability distribution of one overlap except 
for the term "^^a b cQ'abQ^bQ''ca- Putting all terms together we obtain the following expression for the general 
element which involves the probability distribution of three overlaps, 

la^)lai)la^) ^ (^ - '^)(^ -») + »i^ - ^n - l? EabQLEabQUEabQZ .gg^ 

^^^^^^^^ ' (n-2)(n-3)(n-4)(n-5) n{n ~ 1) n{n - 1) n{n - 1) ^ ' 



16 Eab QT^' - 8(n - 2) E. .,. Q'abQTbQi 



n{n - l)(n - 2)(n - 3)(n - 4)(n - 5) 



2(n_8)(n-l) _^Y.abQ'ol' Y.abQZ , Y.abQT' T.abQlb 



(n — 2)(n — 3)(n — 4)(n — 5) n{n — 1) n{n — \) n(n — 1) n{ri — 1) 

n(n — 1) n{n — 1) 
{n - 4)(n - 8) + 8{n - l))(n - l)^ 



(n-2)(n-3)(n-4)(n-5) 
2(n-8)(n-l) 



qkqiQm 



(n-2)(n-3)(n-4)(n-5) 



{qk^m+l + lllm+k + QniQk+l) 



16 8(n-2) 3 

Qm+l+k — 7 7W/ oT? 7T7 ?T *^fe,i ' 



(n - 2)(n - 3)(n - 4)(n - 5) '^™+'+'' („ - 2)(n - 3)(n - 4)(n - 5) 
Where we have introduced the following definitions 



?fe,i ■ 



dqidq2dq3P{qi,q2,q3)q[q2qT and Qt^i^ra,p = dqidq2dq3P{qa,qb,qc,qd)q'aqbQT'1d ■ (^7) 



Finally we get the following expression for the term (p 

{n - 4)(n - 8) + 8(n - l))(n - 1)^ 



[n — 2)(n — 3)(n — 4j(n — 5) 

2(n-8)(n-l) / 2 o N 

(n-2)(n-3)(n-4)(n-5) ^^^ + ^'^^*^ 

16 8(»-2) ^3 



(n-2)(n-3)(n-4)(n-5)'' {n - 2){n - 3){n - A){n - 5) 

For the other term , (q)''', calculations require more effort since one has two deal with terms involving four overlaps, 
but the procedure is similar. We only quote the final result: 

TT4 ^ ((n-l)3(-672 + n(208+(n-24)n)) , 12(n - 1)^(112 + (n - 20)n) , 

^'^' (n-2)(n-3)(n-4)(n-5)(n-6)(n-7) ^^ (n - 2)(n - 3)(n - 4)(n - 5)(n - 6)(n - 7)^^ ^^ ^ ^ 

r^i (-293 + (n - 2) q3^, 1) 



(n - 2)(n - 3)(n - 4)(n - 5)(n - 6)(n - 7) 



12((n-16)(n- 1) 2 20 



92 + 7 ^T7 ?T7 JT7 ET7 ^T7 ^ ^4 



(n - 2)(n - 3)(n - 4)(n - 5)(n - 6)(n - 7) ^ (n - 2)(n - 3)(n - 4)(n - 5)(n - 6)(n - 7) 

8Qli.i-("-3)Qti,i.i) 



4("-2) .,^3 /„ ,^^4 



(n - 2)(n - 3)(n - 4)(n - 5)(n - 6)(n - 7) 

Further calculations lead to a rather intricated expression for Gc in the limit n ^ 0: 

^ 49gf ~ 77g^g2 + gj + 60gig3 - 18g4 + 18QL,i - 24gi Qj,^, - 9 Qj^,^, 

' 133gf - 329g2q2 - 38g| + 240gig3 - 51^4 - 54 Q^ ^^ + 72qiQl^^^^ + 27Qf ^^ ^ 

/,From which the following results for Ac and Be arc obtained: 



_ -49g4 + 77g^g2 - g| - 60gig3 + 18g4 - 18 Qii,i + 24gi Qf^^.i + 9 Qj^^^^^ 
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91gt - 203qh2 + 34g| + SOgigg - 1794 - 18 Qi,i.i + 24gi Q? ^ ^ + 9 q4^^^^^ 

140(92 - qn 

Thus, one has to go to particular structures of the replica matrix Qab to try to obtain a simpler expression for B, A, 
Gc, Be, and Ac. 

The symmetric case 

In the symmetric case, B, the Binder cumulant takes the value 1 in the frozen phase and A vanishes since there 
are no OPF. In all the connected quantities, as weU as in G, both numerator and denominator vanish, thus RE gives 
no information. Despite, we have to stress that, as shown numerically in RS for the spherical SK model, G still takes 
the value 1/3 in the SG phase. 

One-step replica symmetry breaking 

In this particular case the structure of the replica matrix is the following: replicas are distributed into m blocks, so 
that the matrix elements can have two possible values gi if the two replicas belong to the same block and qq if they 
belong to different blocks. Thus, we have for the different terms which appear in the calculation of Gc, A, Ac, B and 
Be-. 

1 >r^ r,k _ (m - l)q'[ + {n - m)q^ 

nin — 1 ^-^ n — 1 

^ ' a,b 

EabcQabQbcQca= (74) 



n(n — l){n — 2) 

91 (m - l)[9i(m - 2) + q^{n - m)] + qa{n - ■m)[2qiqQ{m - I) + q^{n - 2m)] 



n- 1 
1 



n{n- l){n-2) 



j:aAcQabQbcQca= (75) 

qf{m, - l)[qf{m - 2) + ql{n - m)] + 9g(n - m)[2qiqo{m - 1) + q'^{n - 2m)] 

^^1 ' 



- ^ QabQbcQcdQda = (m - l)((m - 3)(rn - 2)qf - 2(3m - 5)(m - n)qlq'^ + 4(2m^ - 3mn + n^)qiql) 



a.b.c.d 



-(3m2 + m(2 - in) + [n - l)2)(m - n)q^ . (76) 



Substituting all these previous expressions into relations K9), ( |66| 



68|) yields, in the limit n — > the expression 



already given in section 111 for the connected quantities Gc (jlS), Ac (|ll) and Be (| 

_ 39-113m + 98TO^ 
" ~ 93-221m + 266m2 ' ^ ' 

In the same way we obtain expressions for Ac and Be which do not depend on the actual values of the elements of 
the replica matrix: 

_ 39-113TO + 98m2 _ 3(31 - 167™ + 182m2) 

^^ - 140m(l - m) ^"""^ ^' - 280m(l - m) ' ^^^> 

For the other disconnected quantities, A and B we obtain an expression which depends not only on m but also on 
the elements of Qafci this is, go ^nd 9i- The expressions are the following: 



(m - 1) m {qf - gg)^ 
3((m- 1) ql -mqlY 



■^ on 1 \ _2 „2^2 ' \ ' 



21 



B 



(m ~ 1) qf — m Qq + 3((to — 1) qf 



mqlf 



2((m -l)q{-m q'^) 



2\2 



(80) 



Comparing A, B to the other connected parameters we find a further dependence on the parameters qa,qi- Compu- 
tations for further steps of RSB are performed in the same way. However for a two-step patter for Qab, we do not 
obtain a simple expression from which it is possible to extract general conclusions. 
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^^ Calculations have been performed under the assumption that there are no singular correlations in the ground state local 

field distribution, which is a direct consequence of dealing with a finite volume system. 
^^ F. Ritort and M. Sales (unpubhshed) 
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